Quantum Hall Effect in Bilayer Graphene: Disorder Effect and Quantum Phase 

Transition 
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We numerically study the quantum Hall effect (QHE) in bilayer graphene based on tight-binding 
model in the presence of disorder. Two distinct QHE regimes are identified in the full energy band 
separated by a critical region with non-quantized Hall Effect. The Hall conductivity around the 
band center (Dirac point) shows an anomalous quantization proportional to the valley degeneracy, 
but the 1/ — plateau is markedly absent, which is in agreement with experimental observation. In 
the presence of disorder, the Hall plateaus can be destroyed through the float-up of extended levels 
toward the band center and higher plateaus disappear first. The central two plateaus around the 
band center are most robust against disorder scattering, which is separated by a small critical region 
in between near the Dirac point. The longitudinal conductance around the Dirac point is shown 
to be nearly a constant in a range of disorder strength, till the last two QHE plateaus completely 
collapse. 

PACS numbers: 73.43.Cd; 72.10.-d: 72.15.Rn 



I. INTRODUCTION 



Since the experimental discovery of an unusual half- 
integer quantum Hall effect (QHE) P, Q in mono- 
layer graphene, the electronic transport properties of 
graphene related materials have been extensively stud- 
ied [a, 4, 5, 6j,0,i,[9iE[ill,[ll,[ll. Recently, bilayer 
graphene is found to show an anomalous behavior in its 
spectral and transport properties, which has attracted 
much experimental and theoretical interest. Theoretical 
studies jil, Q show that interlayer coupling modifies the 
intralayer relativistic spectrum to yield a quasiparticle 
spectrum with a parabolic energy dispersion, which im- 
plies that the quasiparticles in bilayer graphene cannot 
be treated as massless but have a finite mass. Experi- 
ments have shown that bilayer graphene exhibits an un- 
conventional integer QHE [6[. The Landau level (LL) 
quantization results in plateaus of Hall conductivity at 
integer positions proportional to the valley degeneracy, 
but the plateau at zero energy is markedly absent. The 
unconventional QHE behavior derives from the coupling 
between the two graphene layers. The quasiparticles in 
bilayer graphene are chiral and carry a Berry phase 2tt, 
which strongly affects their quantum dynamics. How- 
ever, a detailed theoretical understanding of the uncon- 
ventional properties of the QHE in bilayer graphene tak- 
ing into account of the full band structure and disor- 
der effect is still lacking. As established for a single 
layer graphene [13| and conventional quantum Hall sys- 
tems 



In this work, we carry out a numerical study of the 
QHE in bilayer graphene in the presence of disorder 
based upon a tight-binding model. We reveal that the 
experimentally observed unconventional QHE plateaus 
emerge near the band center, while the conventional QHE 
plateaus appear near the band edges. The unconven- 
tional ones are found to be much more stable to disor- 
der scattering than the conventional ones near the band 
edges. We further investigate the quantum phase transi- 
tion and obtain the phase boundaries Wc for different 
QHE states to insulator transition by calculating the 
Thouless number 15|. Our results show that the un- 



14| . the QHE phase diagram in such a system is 



crucially depending on the topological properties of the 
full energy band, and thus can be naturally determined 
in the band model calculations. 



conventional QHE plateaus can be destroyed at strong 
disorder (or weak magnetic field) through the float-up 
of extended levels toward the band center and higher 
plateaus always disappear first. While the v = ±2 QHE 
states are most stable, the Dirac point at the band center 
separating these two QHE states remains critical with a 
nearly constant longitudinal conductance. 

This paper is organized as follows. In Sec. II, we in- 
troduce the model Hamiltonian. In Sec. Ill, numerical 
results based on exact diagonalization and transport cal- 
culations are presented. The final section contains a sum- 
mary. 



II. THE TIGHT-BINDING MODEL OF BILAYER 
GRAPHENE 

We consider the bilayer graphene composed of two cou- 
pled hexagonal lattice including inequivalent sublattices 
A, B on the bottom layer and A, B on the top layer. 
The two layers are arranged in the AB (Bernal) stack- 
ing 

[3, 

\y\ . as shown in Fig. 1, where B atoms are lo- 
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FIG. 1: (Color online) Schematic of bilayer graphene lattice 
with AB (Bernal) stacking. Bonds in the bottom layer (A, 
B) are indicated by solid lines and in the top layer (^4, B) by 
dash lines. A unit cell contains four atoms: A (white circles), 
B (gray), AB dimer (solid). 



cated directly below A atoms, and A atoms are the cen- 
ters of the hexagons in the^other layer. The unit cell 
contains four atoms A, B, A, and B, and the Brillouin 
zone is identical with that of monolayer graphene. Here, 
the in-plane nearest-neighbor hopping integral between 
A and B atoms or between A and B atoms is denoted 



by JAB 



7 



AB 



— 7o. For the interlayer coupling, we 



take into account the largest hopping integral between B 
and A atoms 7—^ = 71 , and the smaller hopping integral 

between A and B atoms 7^- = 73. The values of these 
hopping integrals are estimated to be 70 = 3.16 eV 
71 = 0.39 eV \M, and 73 = 0.315 eV [20]. 

We assume that each monolayer graphene has totall 



Ly zigzag chains with atomic sites on each chain 



tally 

H 

The size of the sample will be denoted as N = X Ly X 
i^, where 2 is the number of monolayer graphene 

planes along the z direction. In the presence of an applied 
magnetic field perpendicular to the plane of the bilayer 
graphene, the lattice model in real space can be written 
in the tight-binding form: 



iA 



{ij}3 i 

where c|(c|^^), cj(c]g) are creating operators on A and 
B sublattices in the bottom layer, and c|(c^~), cj(c^~) 

iA •' jB 

are creating operators on A and B sublattices in the top 
layer. The sum denotes the intralayer nearest- 

neighbor hopping in both layers, stands for in- 

terlayer hopping between the B sublattice in the bot- 
tom layer and the A sublattice in the top layer, and 
stands for the interlayer hopping between the 

A sublattice in the bottom layer and the B sublattice 
in the top layer, as described above. Wi is a random 
disorder potential uniformly distributed in the interval 
Wi £ [—W/2,W/2]jQ. The magnetic flux per hexagon 



4> = X^o^y ~ ^^^^ ^ an integer. The total flux 
through the sample is N-^, where N = L^Ly/M is taken 
to be an integer. When M is commensurate with or 
Ly, the magnetic periodic boundary conditions are re- 
duced to the ordinary periodic boundary conditions. 



III. RESULTS AND DISCUSSION 

The eigenstates |a) and eigenenergies of the sys- 
tem are obtained through exact diagonalization of the 
Hamiltonian Eq. (1), and the Hall conductivity axy is 
calculated by using the Kubo formula 



'xy 
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a,0 



{a\V^\p){fi\Vy\a)-h.c. 



where S is the area of the sample, Vx and Vy are the 
velocity operators. In Fig. 2a, the Hall conductivity a^y 
and electron density of states are plotted as functions of 
electron Fermi energy Ef for a clean sample {W — 0) at 
system size TV = 96 x 24 x 2 with magnetic flux = |f , 
which illustrates the overall picture of the QHE in the 
full energy band. From the electron density of states, 
we can see the discrete LLs. We will call central LL at 
Ef = the n = LL, the one just above (below) it 
the n = 1 {n — —1) LL, and so on. According to the 
behavior of axy, the energy band is naturally divided 

into three different regimes. Around the band center, 

2 

the Hall conductivity is quantized as axy = ^^X' "^here 
v = kgs with k an integer a nd p s = 2 for each LL due to 
double- valley degeneracy 0, [l3| (the spin degeneracy will 
contribute an additional factor 2, which is omitted here). 

With each additional LL being occupied, the total Hall 

2 

conductivity is increased by 5s 77- This is an invariant 
as long as the states between the n-th and {n — l)-th 
LL are localized, axy = at the particle-hole symmetric 
point Ef = 0, which corresponds to the half-filling of 
the central LL. However, there is no axy — quantized 
Hall plateau. These anomalously quantized Hall plateaus 
agree with the results observed experimentally in bilayer 
graphene 0]. 

The Hall conductivity near the band edges, however, is 
quantized as axy = with k an integer, as in the con- 
ventional QHE systems. Remarkably, around Ef — ±70 
(within a narrow energy region Ai? ~ Q.470), there are 
two critical regions which separate the unconventional 
and conventional QHE states, where the Hall conduc- 
tance quantization is lost. These crossover regions also 
correspond to a novel transport regime, where the Hall 
resistance changes sign and the longitudinal conductiv- 
ity exhibits metallic behavior. The singular behavior of 
the Hall conductivity in the crossover regions is likely to 
originate from the Van Hove singularity in the electron 
density of states at i? = limit. In Fig. 2b, the quantiza- 
tion rule of the Hall conductivity in this unconventional 
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FIG. 2: (Color online) (a) Calculated Hall conductivity and 
electron density of states in the full energy band for mag- 
netic flux = ll or M=48, and (b) the Hall conductivity 
near the band center for </> = , |f- and The disorder 
strength is set to W = and A'^ = 96 x 24 x 2 in all cases. 
Inset: Hall conductivity at the band center. Here, the spin 
degrees of freedom are omitted, so (7s = 2 and Qs — 1 for the 
unconventional and conventional regions, respectively. 
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FIG. 3: (Color online) Unconventional Hall conductivity as a 
function of electron Fermi energy near the band center for four 
different disorder strengths each averaged over 400 disorder 
configurations. Inset: conventional Hall conductivity near 



the lower band edge. 
iV = 96 X 24 X 2. 



Here, 



and the sample size is 



region for three different strengths of magnetic flux is 
shown. With decreasing magnetic flux from (f> — ^ to 
ll, more quantized Hall plateaus emerge following the 
same quantization rule as the gap between the LLs is 
reduced. 

Now we study the effect of random disorder on the 
unconventional QHE in bilayer graphene. In Fig. 3, the 
Hall conductivity around the band center is shown as 
a function of Ef for four different disorder strengths at 



FIG. 4: (Color online) Calculated Hall conductivity with 
weaker magnetic flux = ||-, and ^ for four different 
disorder strengths each averaged over 400 disorder configura- 
tions. Here, the sample size is iV = 96 x 24 x 2. 



system size = 96 x 24 x 2 with magnetic flux cj) — 
We can see that the plateaus with v = ±10, ±6 and ±2 
remain well quantized at W = 0.5. We mention that the 
h' = ±4, ±8 plateaus are unclear at this relatively weak 
disorder strength because of very small plateau widths 
and relatively large localization lengths (the critical Wc 
for each plateau will be obtained based on our larger size 
calculations of the Thouless number as presented later). 
With increasing W, higher Hall plateaus (with larger \iy\) 
are destroyed first. At W = 2.0, only the = ±2 QHE 
remain robust. The last two plateaus = ±2 eventually 
disappear around W ^ 3.2. For comparison, the QHE 
near the lower band edge is shown in the inset, where all 
plateaus disappear at a much weaker disorder strength 
W > 1.0. This clearly indicates that under the same con- 
ditions, the unconventional QHE around the band center 
is much more stable than the conventional QHE near the 
band edges. Clearly, after the destruction of the con- 
ventional QHE states near the band edge, these states 
become localized. Then the topological Chern numbers 
initially carried by these states will move towards band 
center in a similar manner to the single-layer graphene 
case 



13|. Thus we observe that the destruction of the 



unconventional QHE states near the band center is due 
to the float-up of extended levels. 

To study the fate of the IQHE at weak magnetic 
field limit, we reduce the strength of magnetic field. In 
Fig. 4, the Hall conductivities around the band center 
with weaker magnetic fiux (j) — and are 

shown for different disorder strengths and system size 

= 96 x 24 X 2. In Fig. 4a, a lot more well quantized 
Hall plateaus emerge for a clean sample(iy — 0), if we 
compare them with the results in Fig. 2b. In Fig. 4b, 4c 
and 4d, we can see that with the increasing of the dis- 
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order strength W, Hall plateaus are destroyed faster for 
the system with weaker magnetic flux 0. At = 2.0, 
the most robust Hall plateaus at i' ~ ±2 remain well 
quantized for magnetic flux (Z* = |f and j^, however, 
they already disappear for weaker magnetic flux </) = ^ . 
Our flux 2tt/M in each hexagon the magnetic field is 
S ^ 1.3 X lO^/M TeslapH. Thus the weakest B we 
used is about 451 Tesla. This is a very large magnetic 
field comparing to the experimental ones around B ^ 40 
Tesla. However, the topology of the QHE and how they 
disappear with the increase of the disorder strength W 
remain to be the same as the stronger B cases as demon- 
strated in Fig. 4a-4d. Thus, we establish that the ob- 
tained behavior of QHE for bilayer graphene will survive 
at weak B limit. 

We further investigate the quantum phase transition 
of the bilayer graphene electron system. In order to de- 
termine the critical disorder strength Wc for the different 
QHE states, the Thouless number g is calculated by using 
the following formula [3], 
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AE 

dE/dN 



Here, AE is the geometric mean of the shift in the en- 
ergy levels of the system caused by replacing periodic 
by antiperiodic boundary conditions, and dE/dN is the 
mean spacing of the energy levels. The Thouless number 
g is proportional to the longitudinal conductance G. In 
Fig. 5, we show some examples of calculated Thouless 
number for a relatively weak flux 4> — ^ and some dif- 
ferent disorder strengths to explain how quantum phase 
transitions and the related phase boundaries Wc are de- 
termined. In Fig. 5a, the calculated Thouless number g 
and Hall conductivity a^y as a function of i?j at a weak 
disorder strength W = 0.2 are plotted. Clearly, each 
valley in Thouless number corresponds to a Hall plateau 
and each peak corresponds to a critical point between 
two neighboring Hall plateaus. We can also call the first 
valley just above (below) Ef = the = — 2 (i^ = 2) 
QHE state, the second one the v — —4 {v — 4) state, 
and so on, as same as the Hall plateaus. In Fig. 5b-5d, 
we see that with increasing W, higher QHE states (val- 
leys) are destroyed first. AtW = Wc = 1.0 (see Fig. 5b), 
the valleys with v — ±12 disappear, which correspond 
to the destruction of the — ±12 Hall plateau states. 
Therefore, Wc = 1.0 is the critical disorder strength, at 
which the v = ±12 plateau states change to an insulat- 
ing phase. At W = Wc = 1-3 (see Fig. 5c), the valleys 
with V = ±8 disappear, which indicates the destruction 
of the V — ±8 QHE states and their transition into the 
insulating phase. When W — Wc — 3.2 (see Fig. 5d), the 
most stable QHE states with v = ±2 eventually disap- 
pear, which indicates all QHE phases are destroyed by 
disorder. All the phase boundaries Wc between the dif- 
ferent QHE states are determined in the same manner 
and tabulated in Table 1. 
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FIG. 5: (Color online) (a)-(c) Calculated Thouless number 
and Hall conductivity for three different disorder strengths, 
and (d) Thouless number for other three disorder strengths, 
each data point being averaged over 400 disorder configura- 
tions. Here, and the sample size are taken to be 
TV = 96 X 48 X 2 and ./V = 96 x 24 x 2 in the calculations of 
Thouless number and Hall conductivity, respectively. 



We now focus on the region around Ef = 0. In Fig. 6, 
we show the Thouless number for some different disorder 
strengths at system size N — 96 x 24 x 2 and magnetic 
flux (f> — j^- We can see that the Thouless number shows 
a central peak at Ef = 0. With increasing the disorder 
strength, the width of the peak increases and its height 
remains nearly unchanged. This behavior may suggest 
an interesting effect that the extended states originally 
sited at the critical point Ef = splits in the presence 
of disorder. However, the splitting is too small to in- 
duce two separated peaks in the Thouless number for the 
present sample sizes we can approach. Instead, it leads 
to a widened peak of unreduced height. This behavior 
also indicates that the critical longitudinal conductance 
in a small finite region near Ef = is almost constant 
about 2e^/h according to the proportionality of Thou- 
less number to longitudinal conductance. We have also 
confirmed this conclusion by direct Kubo formula calcu- 
lation, in which the system size that can be approached 
is however much smaller. 



Hall plateaus index critical point Wc 



V = ±12 


1.0 


V = ±10 


1.2±0.1 


i/ = ±8 


1.3±0.1 


= ±6 


1.6±0.1 


= ±4 


1.7±0.1 


= ±2 


3.2 



TABLE I: The phase boundaries Wc. for the different Hall 
plateaus. 
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thus one will not be able to see the important physics of 
opposite Chern numbers annihilating each other to de- 
stroy the IQHE[ii. 
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FIG. 6: (Color online) (a) Thouless number for five different 
disorder strength, each point being averaged over 400 disorder 



configurations. 
96 X 48 X 2. 



Here, 



7?- and the sample size is A'' 



IV. SUMMARY 

In summary, we have numerically investigated the 
QHE in bilayer graphene based on tight-binding model 
in the presence of disorder. The experimentally ob- 
served unconventional QHE is reproduced near the band 
center. The unconventional QHE plateaus around the 
band center are found to be much more stable than the 
conventional ones near the band edges. Our results of 
quantum phase transition indicate that with increasing 
disorder strength, the Hall plateaus can be destroyed 
through the float-up of extended levels toward the band 
center and higher plateaus always disappear first. At 
W = Wc = 3.2, the most stable QHE states with v = ±2 
eventually disappear, which indicates transition of all 
QHE phases into the insulating phase. A small crit- 
ical region is observed between the v — ±2 plateaus, 
where the longitudinal conductance remains almost con- 
stant about 2e^//i in the presence of moderate disorder, 
possibly due to the splitting of the critical point origi- 
nally sited at Ef = 0. We mention that in our numerical 
calculations, the magnetic field is much stronger than 
the ones one can realize in the experimental situation, as 
limited by current computational ability. However, the 
phase diagram we obtained is robust and applicable to 
weak field limit since it is determined by the topologi- 
cal property of the energy band as clearly established for 
single layer graphene [l3| and conventional quantum Hall 
systems IJ]. We further point out that the continuum 
model can not be used to address the fate of the quan- 
tum Hall effect in strong disorder or weak magnetic field 
limit. Because in such a model, both the band bottom 
and band edge are pushed to infinite energy limit, and 
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